Is timing noise important in the gravitational wave detection of neutron stars? 
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In this paper we ask whether the phenomenon of timing noise long known in electromagnetic 
pulsar astronomy is likely to be important in gravitational wave (GW) observations of spinning- 
down neutron stars. We find that timing noise is strong enough to be of importance only in the 
young pulsars, which must have larger triaxialities than theory predicts for their GW emission to be 
detectable. However, assuming that their GW emission is detectable, we list the pulsars for which 
timing noise is important, either because it is strong enough that its neglect by the observer would 
render the source undetectable, or else because it is a measurable feature of the GW signal. We also 
find that timing noise places a limit on the observation duration of a coherent blind GW search, and 
suggest that hierarchical search techniques might be able to cope with this problem. Demonstration 
of the presence or absence of timing noise in the GW channel would give a new probe of neutron 
star physics. 



PACS numbers: 04.30.Db, 04.40.Dg, 04.80.Nn, 97.60. Gb 



I. INTRODUCTION 

Spinning triaxial neutron stars may provide a source 
of detectable gravitational radiation for the new genera- 
tion of gravitational wave interferometers |l| . The known 
pulsar population will be targeted in forthcoming GW 
searches, as their locations are known to high accuracy 
and their rotational properties have been studied in some 
detail, greatly aiding the search. The GW amplitudes 
of these sources are likely to be low, so that only by 
coherently accumulating signal over long time intervals 
(months to years) is there any hope of making a posi- 
tive detection [2| • This accumulation can be achieved by 
the process of matched filtering, where the noisy detec- 
tor output is multiplied by a template waveform which 
remains in phase with the GW signal to better than one 
radian over the entire observation span, and the resulting 
product integrated. 

The spin frequencies of almost all pulsars are observed 
to gradually decrease, presumably because of the loss of 
energy caused by electromagnetic and (hopefully) GW 
emission [3j- Apart from occasional glitches, this fre- 
quency change is gradual, so that pulsar physicists can 
model the rotation phase using only a few terms of a 
Taylor series of the form 



$ = 2tt / f + f t + fot 2 /2 + ...dt. 



(1) 



However, accurate radio timing observations reveal a 
small irregularity in this spin-down known as timing 
noise, occurring on long timescales, comparable with the 
likely GW observation timescales Q . This noisy behavior 
of the electromagnetic signal may well be present in the 
GW signal too, as almost all parts of a neutron star are 
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believed to be coupled together strongly on a timescale 
of seconds or less [5J . 

If present in the GW signal the timing noise would 
result in a phase difference between the idealized Taylor 
expansion and the real noisy GW signal which grows in 
time (see Equation below). It follows that there would 
exist a timescale Tdccoherencc at which a simple Taylor 
series of the above form will drift out of phase with the 
template by one radian, leading to a complete loss of 
signal-to-noise and preventing detection. This could be 
prevented by the GW observer going to the trouble of 
using the electromagnetically detected phase to generate 
the GW template. A method for achieving this has been 
devised by Pitkin & Woan for the Crab pulsar @ ■ 

However, the timing noise should not be regarded sim- 
ply as a nuisance, complicating the detection process. As 
described in section [HI the strength of the timing noise 
can be used to probe stellar structure, giving a unique 
insight into neutron star physics. The accuracy with 
which the phase of the GW signal can be extracted from 
the noisy data stream is proportional to one over the 
signal-to-noise ratio |7J, and so the GW phase measure- 
ment error decreases as T -1 / 2 , where T is the duration 
of the coherent observation. It follows that there exists a 
timescale ^detectable at which the strength of the timing 
noise in the GW data stream can be measured. 

With these remarks in mind, we will pose and answer 
the following three questions in this paper: 

• For each known pulsar, on what timescale 
^decoherence would the timing noise, if present in its 
GW signal, cause a simple Taylor series search tem- 
plate to drift out of phase with the actual signal, 
leading to a complete loss of signal-to-noise? 

• For each known pulsar, on what timescale 
^detectable would the strength of GW timing noise 
be measurable? 

• Is timing noise likely to be important when per- 
forming blind searches, i.e. GW searches for neu- 



tron stars not currently observed as pulsars? 



tween the GW and EW phase residuals 



In all cases we have chosen to restrict our attention to 
observation durations of three years or less. This is partly 
because the computational costs of long coherent obser- 
vations are very large 8] , and partly because the interfer- 
ometers are likely to be upgraded significantly over such 
a timescale, in which case analysis of a smaller set of 
newer, higher quality data might be more profitable. 

The structure of this paper is as follows. We begin 
by summarizing the physical significance of timing noise 
in section [H] In section IIIII we assemble formulae nec- 
essary to estimate T dccoh oionco and T dGt octabio- In section 
II VI we describe how we estimate the strength of timing 
noise in the known pulsar population. In section we 
provide both upper bounds and somewhat more realis- 
tic estimates of the GW amplitudes of these pulsars. In 
section IVll we estimate the decoherence timescales of the 
pulsars. In section IVlII we estimate possible timescales 
for detectability of timing noise in the GW signal. In 
section I Villi we briefly discuss the importance of timing 
noise for blind searches. We summarize our results in 
section IIXI 



A$ 



GW 



low 



A$EW- 



(3) 



II. THE PHYSICAL SIGNIFICANCE OF 
TIMING NOISE 

Theoretically, the origin of timing noise is not under- 
stood. Eight different models of timing noise were consid- 
ered in J3] , where the statistical properties of the timing 
residuals of each model were calculated and compared 
with electromagnetic observations. Some of the models 
performed better than others, but a definitive identifica- 
tion of the basic mechanism at work proved impossible. 

Nevertheless, it is instructive to consider three rather 
general scenarios, without linking them to any particular 
model of timing noise. In the first the whole star, in- 
cluding the pulsation-producing magnetosphere, rotates 
as a rigid body, so that the relative phase of the GW 
and electromagnetic wave (EW) is constant in time. In 
the second, the timing residuals represent a purely mag- 
netospheric phenomenon, where the location of pulsa- 
tion production, at a height of several neutron star radii 
above the surface, wanders randomly in longitude, with- 
out there being any corresponding variation in the ro- 
tational phase of the star. In the third scenario, the 
variation in phase is ascribed to a weak random angular 
momentum exchange between the part of the star tied to 
the EW emission and the part tied to the GW emission. 
Denoting the moment of inertias of these two parts by 
Jew and /qw, and the departure of their angular veloci- 
ties from a smooth spindown law by Afisw and A£!gw> 
conservation of angular momentum then demands 

IewASIew + JgwA^gw = 0. (2) 

Integrating twice and rearranging gives a relationship be- 



(As we are considering triaxial neutron stars, the GW 
emission is at twice the rotation frequency |l0| . so that 
the quantity A$qw in Equation © is strictly only one- 
half of the GW phase residual). In this scenario, the 
anti- correlation of the gravitational and electromagnetic 
phase residuals would then allow us to probe the rela- 
tive moments of inertia of the two parts of the star, as 
well as telling us that they are loosely coupled on the 
timescale of the phase wandering. For instance, if these 
two parts were the crust and fluid core, we would expect 
Iew/Igw-10- 2 O. 

Theoretically, all parts of the star (apart from any su- 
perfluid pinned to the crust) are expected to be coupled 
strongly on timescales of the order of seconds y| , orders 
of magnitudes smaller than the timescales of years on 
which the phase residuals vary. It therefore seems most 
likely that the first of the above three scenarios is cor- 
rect. However, given the potential importance of this 
phase information, the following strategy suggests itself: 
Wherever possible, the GW analysis of a known pulsar 
should use the observed pulsar phase residuals A$em 
to demodulate the GW detector output prior to using a 
smooth spindown template of the form of equation Q. 
A way in which such an analysis could be performed for 
the Crab pulsar has been discussed recently by Pitkin 
&; Woan |g. The demodulation itself should correct for 
GW phase wandering of the form 



A$gw = aA$ 



EW, 



(4) 



with a range of different a values used. The three scenar- 
ios above correspond to the cases a = 1,0, —Iew/Igw, 
respectively. Measurement of a value of a other than 
unity would challenge the standard picture of neutron 
star structure. 

Of course, the addition of an extra parameter to an 
already computationally expensive search should only be 
done if the extra parameter is likely to have a significant 
affect upon the signal — either by making what would oth- 
erwise be an undetectable signal detectable, or by signif- 
icantly improving the signal-to-noise of an already de- 
tectable one. It is precisely these issues that we address 
in this paper, to assess whether or not GW observers 
need go to the trouble of allowing for timing noise in 
their analyzes. 



III. METHOD OF CALCULATION 

Timing noise can be characterized as a random walk 
in one or a combination of the rotation phase, frequency 
or spindown rate Q- These three idealized behaviors 
are known as phase noise, frequency noise and spindown 
noise. Pulsar physicists quantify the strength of timing 



noise by < A$| w > 1 / 2 , the root mean square phase resid- 
ual separating the actual signal and the best Taylor series 
approximant to it, typically containing terms up to and 
including the first frequency derivative. For the three 
idealized forms of random walk this rms phase residual 
grows as: 



<a$| w > 1 / 2 =a : t"/ 2 , 



(5) 



where n = 1,3,5 for phase noise, frequency noise and 
slow-down noise, respectively. The corresponding GW 
timing noise in our model is then: 



<A$ 2 , w > 1 / 2 =| a |fcT"/ 2 . 



(6) 



A. Calculation of the decoherence timescale 

Taking our criterion for total decoherence to be 
<A$ 2 3W > 1 / 2 = 1 radian, Equation (JBJ gives 



' decoherence 



Ijn 



(7) 



B. Calculation of the detection timescale 

A constant frequency periodic GW source is described 
by 7 parameters: two angles giving its direction in the 
sky, its frequency, the phase at time t — 0, two angles 
specifying the orientation of its spin axis, and its ampli- 
tude |l2j. When searching for GWs from a known pulsar, 
the first four of these are known to high accuracy, so only 
the last three need be fit for. 

When timing noise is present in the GW signal (at 
least) two possible strategies present themselves. Firstly, 
we could imagine taking the total data stream and break- 
ing it up into a small number of equal duration blocks, 
sufficiently short that each block can be analyzed without 
worrying about the effects of timing noise. The timing 
noise would then manifest itself as a smooth variation in 
the phase at t = from block to block. This method has 
the advantage that we need make no assumptions about 
the relative phase of the EW and GW timing noises, i.e. 
does not assume a correlation of the form of Equation 0} , 
but the disadvantage that if the number of data blocks 
is large the signal-to-noise ratio for each will be small, 
leading to inaccurate phase measurements. 

We will pursue a second, somewhat simpler strategy in 
this paper. We will assume that the EW and GW timing 
noises are correlated as in Equation Q . Then a must be 
added as an eighth search parameter. This method has 
the advantage that the entire signal-to- noise of the obser- 
vation can be brought to bear on the analysis. Of course, 
if the EW and GW timing noises are not correlated we 
would not expect any value of a to significantly maximize 
the signal-to-noise ratio, indicating that we must use the 
first method described above. 



We present below a simple estimate of how long an 
observation time is required to measure a, i.e. to obtain 
a value significantly different from zero. 

The error in measuring the phase of a periodic GW 
signal of signal-to-noise ratio p is given approximately by 



<5$ 



GW 



N 
P 



(8) 



where N is a small number which depends upon the num- 
ber of parameters to be fit for yfl. The signal-to-noise 
grows as T 1 ' 2 ; we can make this time dependence ex- 
plicit by writing 



8$ 



N 



GW 



Pi 



T l/2' 
-yr-fyrs 



(9) 



where pi_ yr is the signal-to-noise that would be attained 
for a one year observation, and T yrs is the total observa- 
tion duration in years. 

We will take as our criterion for the detectability of 
timing noise in the GW signal that the rms phase wan- 
dering is equal to the GW phase measurement error, i.e. 
< A$q W > 1 / 2 w 5$gw- Combining Equations 10 and 
© gives the time T equa i at which this occurs 



' equal, yrs 



N 



\a\kpi 



2/(n+l) 



(10) 



Of course, it is also necessary that the GW signal be 
detectable. Let us set the minimum signal-to-noise ratio 
for detection to some value /9 m i n - Define this signal-to- 
noise to be attained for an observation duration Tqw- 
Using the result that the signal-to-noise grows as T 1 / 2 
we have: 



?bw = 



Pl-yr 



(11) 



The timing noise is only detectable when its magnitude 
exceeds the GW phase measurement error (i.e. T > 
Tequai) and the GW signal is detectable (i.e. T > Tqw)- 
It therefore follows that the presence of timing noise can 
be detected in the GW data stream after an observation 
time Tdetectabie equal to the maximum of these two val- 
ues: 



T, 



detectable 



= max(T cqua i, T G w)- 



(12) 



We would like to apply the above formulae to the entire 
pulsar population, estimating T dot ectabie and T dcconcrcnco 
for each known pulsar. In order to make use of the for- 
mulae, we need to obtain estimates of the timing noise, 
parameterized by k and n, and also of the GW signal 
strength, parameterized by pi- yr . 



IV. THE TIMING NOISE 

Ideally we would use observationally derived timing 
noise parameters k and n for the entire pulsar popula- 
tion. Unfortunately, a literature search shows that only 



TABLE I: Papers containing radio pulsar timing noise data. 
The number of pulsars provided with /gw > 10 Hz or greater 
is indicated. 



Reference 


Number of pulsars 


Cordes & Downs (1985), [171 


1 


Cordes & Helfand (1980), [l£] 


2 


D'Alessandro et al. (1995), [191 


2 


Cordes et al. (1988), [201 


1 


Manchester et al. (2001), 1151 


10 


Morrisat et al. (2002), [161 


14 


Kramer et al. (2003), [14] 


23 



a minority of pulsars have been timed in sufficient de- 
tail to allow for estimation of the timing noise behavior. 
However, on the basis of the available data, Dewey & 
Cordes [l3| have obtained a fitting formula that allows 
estimation of the strength of the timing noise of a pulsar 
in terms of its period P and period derivative P. We 
will therefore pursue the following strategy: For pulsars 
whose timing noise strength has been measured, we will 
use the observational data in our analysis. For pulsars 
whose timing noise hasn't been measured (or at least 
weren't provided by our literature search), we'll use the 
fitting formula of Dewey & Cordes. 



A. Pulsars of measured timing noise 



A number of authors have published rms phase resid- 
uals < A$| w >U 2 for pulsars, typically after having fit 
the timing data to a Taylor expansion including the fre- 
quency and its first time derivative, over an interval of 
several years. The parameter k can then be derived us- 
ing the published values of < A$| w > 1 / 2 , T and n in 
Equation ©: 



_ <A$| W >V 2 

ft rpn/2 



(13) 



Unfortunately, in only a small subset of these observa- 
tions was it possible to identify the nature of the walk, 
i.e. fix the value of n. For the sake of defmiteness, we 
will therefore present results for n — 3; we have repeated 
calculations for n = 1 and n = 5 and found results which 
are not very different, with the n = 3 results being inter- 
mediate between the two other sets. 

The results of out literature search for timing noise 
data are sunrmarized in Table [I] For reasons of con- 
venience, we only selected pulsars of spin frequencies 
greater than 5 Hz, corresponding to GW frequencies 
greater than 10 Hz, as stars spinning more slowly than 
this will certainly not be of GW interest. The references 
of table U then provided a total of 53 pulsars. Note that 
the majority of the pulsars selected come from the Parkes 
Multibeam Study p^ll5lllq |. which was specifically de- 
signed to find young pulsars. 



B. Pulsars whose timing noise hasn't been 
measured 



In producing their fitting formula for timing noise, 
Dewey & Cordes make use of the 'activity parameter', 
defined as the logarithm of the rms phase residual for 
the pulsar in question divided by that of the Crab: 



A = lop 



<A$| W > 1 / 2 

< A *l W ,Crab> 1/2 



(14) 



The fitting formula is 

J 4 = -1.41ogP + 0.81ogP_i5-3.31. (15) 

where P is the pulsar's spin period in seconds and P 
is the dimensionless period derivative divided by 1CP 15 . 
The Crab is observed to display frequency-type timing 
noise of the form: 



<A*BW,Crab> 



1/2. 



0.24 radians T y 3 / S 2 



(16) 



(see |4|). It follows that the activity parameter only 
makes sense for (Crab-like) frequency-type noise, as A 
would otherwise be a function of time. We will there- 
fore assume that all pulsars of measured timing noise dis- 
play frequency-type noise. The above three equations can 
then be combined to give an estimate of the rms timing 
noise in any pulsar of known period and period deriva- 



tive. Using <A$ 2 =W > 1/2 / < A$ 



EW, Crab 



> 



1/2- 



k/k 



in Equation 114fl then allows calculation of k: 



k = fccrablO- 4 . 



Crab 



(17) 



Application of this formula to pulsars of known timing 
noise shows that there exists a scatter of about an order 
of magnitude in the measured activity parameter values 
about the predicted ones [J • This uncertainty should be 
borne in mind when reading the results sections of this 
paper — individual pulsars may be more or less noisy than 
assumed. 

We used the above prescription for the 1182 pul- 
sars listed in the Australia Telescope National Facility 
database [21j whose periods, period derivatives and dis- 
tances from Earth were known, excluding the 53 of mea- 
sured timing noise described in section IIV Al i.e. the 
fitting formula was used for 1182 — 53 = 1129 pulsars. 



V. GRAVITATIONAL WAVE AMPLITUDES 

We also need to estimate the GW amplitudes of the 
pulsar population, giving the result in the form of the 
signal-to-noise for a one year coherent integration. The 
GW amplitude of a triaxial star a distance r from Earth 
is given by: 



"='15 



1/2 



Gsn 2 



A/, 



(18) 



where the normalization comes from averaging over all 
possible spin orientations of the source and A/ is the 
difference between the two principal moments of inertia 
in the plane orthogonal to the spin 2]. This is more 
conveniently expressed as a dimensionless number 



A£ 

T' 



(19) 



where / is the moment of inertia of the non-rotating star; 
we will refer to e as the ellipticity. Of course, the wave 
amplitudes at Earth are unknown as e is unknown. An 
upper bound is often placed on these quantities by as- 
suming that all the the spindown kinetic energy lost by 
the star is converted into GW energy, at fixed moment 
of inertia. In general the GW luminosity is given by 



E = 



32G 

5c5" 



n 6 (Ai) 2 



(20) 



while the kinetic energy is 7T2 2 /2. Combining these re- 
sults gives an upper bound on h in terms of the pulsar's 
spin period, period derivative and distance from Earth: 



2 GIP 



1/2 



C*P 



(21) 
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FIG. 1: Upper bounds on GW amplitudes assuming 100% 
conversion of spindown energy into GWs. The noisecurves 
are for a one year observation, and were produced using the 
fitting formulae of |22| . The ten red circles indicate those 
pulsars which have Tdecoherenca < 10 years; see Section IVT1 



with a corresponding ellipticity 

5c 5 PP 3 



^spindown 



32(27r) 4 GI 



1/2 



(22) 



The results of such a calculation for all pulsars of known 
period, period derivative and distance from Earth are 
shown in Figure ^ 

The ellipticities required to produce this level of spin- 
down are of order 10 -8 for the millisecond pulsars. De- 
tailed modeling suggests that such an ellipticity could 
well be produced by strains in the neutron star crust 
|23|. However, for the younger pulsars, which include all 
those which lie above the noisecurves of the first genera- 
tion detectors, the required ellipticities are much larger, 
around 10~ 4 . Such a value is almost certainly unphysi- 
cally large, suggesting that GW play little role in the en- 
ergy budget of young neutron stars. With this in mind, 
in Figure |21 we show another set of estimated wave am- 
plitudes, repeating the above calculation, but imposing a 
cut-off in e of 10~ 7 , i.e. setting e = max(e sp i n down, 10 -7 ). 
The value e = 10~ 7 was chosen as it is at the upper- 
end of what is considered plausible on physical grounds 
J23J. The wave amplitudes calculated for the millisec- 
ond pulsars are mostly unchanged, but the amplitudes 
for the young pulsars have been reduced by several or- 
ders of magnitude (most have fallen off the bottom of 
the figure). In the following sections we will present re- 
sults assuming ellipticities ranging from e S pindown down 
to 1CP 7 . (It is probable that in reality pulsar ellipticities 
are even smaller than 10~ 7 , but, as we shall see, we need 
not consider smaller values for the purposes of examining 
timing noise). 
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FIG. 2: Upper bounds on GW amplitudes imposing at cut-off 
of 10 -7 in e. 



VI. RESULTS: THE DECOHERENCE 
TIMESCALE 



We will now insert the timing noise strengths calcu- 
lated in section lTvl into Equation|7|to obtain estimates of 
the timescale on which timing noise would cause a simple 
Taylor series to completely decohere from the GW signal, 
assuming that the latter is perfectly locked in phase with 
the EW signal (i.e. a = 1 in Equation 0]). These esti- 
mates are completely independent of the GW strength — 
they depend only upon the strength of the timing noise. 

We find that 10 pulsars have decoherence timescales 



TABLE II: Pulsars with Tdecoherence < 3 years and Tgw < 
3 years, assuming the Initial LIGO noisecurve, 100% conver- 
sion of spindown energy to GW energy, and a — 1. 



Name 


-*- dccohcrcncc/ yCoXS 


Tgw/ years 


J0205+6449 

J0537-6910 

B0531+21 


2.0e+00 
l.le+00 
2.6c+00 


2.8e+00 
2.3e+00 
1.3e-02 



TABLE III: Pulsars with Tdecoherence < 3 years and Tgw < 
3 years, assuming the Advanced LIGO noisecurve, 100% con- 
version of spindown energy to GW energy, and a — 1 



Name 


J- dccohcrcncc/ y63XS 


Tgw/ years 


J0205+6449 


2.0e+00 


1.4e-04 


B0531+21 


2.6e+00 


1.9e-06 


J0537-6910 


l.le+00 


6.2e-03 


B0540-69 


9.5e-01 


4.1e-03 


J1124-5916 


1.9e+00 


3.7e-03 


B1509-58 


1.4e+00 


2.0e-03 


J1617-5055 


2.5c+00 


l.le-03 


J1930+1852 


1.9e+00 


3.3e-03 


J2229+6114 


2.6e+00 


1.7e-03 



of less than three years. To gain some insight, we show 
their positions in the GW amplitude-frequency plot in 
Figure ^ as open circles, assuming 100% conversion of 
spindown energy into GWs. As expected, most of these 
are young pulsars — only one has a GW frequency of less 
than 10 Hz. None of them are millisecond pulsars; the 
much weaker timing noise of the latter class of star leads 
to their having decoherence timescales many orders of 
magnitude longer that those of the young pulsars. 

However, these estimates are only of interest if the GW 
signal is itself detectable. To this end, we will present re- 
sults first assuming 100% conversion of spindown energy 
into GWs, and then assuming various bounds upon the 
allowed ellipticity. 



for which timing noise could prevent detection. In fact, 
detection will be prevented only if Tdecoherence < Tgw- 
For the 100% conversion of spindown energy to GW en- 
ergy assumed in Tables |H] and Hill this is the case only 
for the Initial LIGO observation of pulsars J0205+6449 
and J0537-6910. 

Note that only one pulsar has a predicted decoherence 
timescale of less than one year, namely pulsar B0540- 
69, which has Tdecoherence = 0.95years. We therefore 
see that, regardless of the GW signal strength, timing 
noise can only ruin a Taylor-series based GW detection 
for timescales of order a year or longer. 

Motivated by the discussion of section[Hj if we set |a| = 
10~ 2 in Equation J7J) then, for frequency type noise, the 
decoherence timescales will be greater than those of Table 
IIIII by a factor of 100 2 / 3 ~ 21.5; in this case there are no 
known pulsars with Tiecohercncc < 3 years. 



B. Decoherence timescale assuming various upper 
bounds on the triaxiality 

Putting an upper bound on the ellipticity will have 
no effect on the decoherence timescales but will increase 
the GW detection timescales, so that for each pulsar in 
the Tables ITT1 and IIIII for some assumed e max , the con- 
dition Tdecoherence < Tqw would be satisfied and tim- 
ing noise would prevent GW detection. For instance, 
putting e max = 6x 10~ 7 the Advanced LIGO observation 
of the Crab satisfies the condition. However, repeating 
the analysis for the more physically plausible case where 
£max = 10~ 7 , we find that the GW detection timescales 
of all pulsars in the table exceed 3 years, even for the 
Advanced LIGO noisecurve. 

To sum up, timing noise could render as many as 3 
pulsars undetectable using simple Taylor series for Initial 
LIGO, and as many as 9 for Advanced LIGO, but GW 
emission from these pulsars is only detectable if they have 
ellipticities in excess of those expected theoretically. 



A. Decoherence timescale assuming 100% 

conversion of spindown energy into GWs; probably 

unphysical 

As stated above, we find that 10 pulsars have decoher- 
ence timescales of less than three years. However, even 
when making the optimistic assumption that 100% of 
the spindown energy is being converted into GW energy, 
some of these pulsars have GW detection timescales of 
more than three years, and are therefore of no interest. 
To allow for this, the names, Tdecoherence and Tqw values 
for the pulsars which satisfy Tdecoherence < 3 years and 
Tgw < 3 years are given in Table [H] assuming the Ini- 
tial LIGO noisecurve, and the corresponding results for 
Advanced LIGO are given in Table llTTl 

Tables [H] and IIIII represent the maximal set of pulsars 



VII. RESULTS: THE DETECTION TIMESCALE 

We will now combine the timing noise estimates of sec- 
tion IIVI with the GW amplitudes of section to calcu- 
late the Tietectabic values for all known pulsars. We will 
present results assuming both the Initial and the Ad- 
vanced LIGO noisecurves. In the results that follow we 
have set the parameter N = 3 in Equation [§| Cutler 
et al. (2003) [Tfl have shown that for LISA observations 
at least, N < 1.5 for 40% of the possible spin orienta- 
tions of the source, so our setting N to twice this value is 
surely a conservative estimate of the phase error. Also, 
we have set p mm — 5 in Equation Ijllj l: for the directed 
searches considered here this will give an acceptably low 
false alarm rate (Reference Q, Eqn. 1.4). 

Again, we will divide our results into two parts. In the 
first part we will assume complete conversion of spindown 



TABLE IV: Pulsars with Tdctoctabio < 3 years, assuming 100% 
conversion of spindown energy into GW energy, the Initial 
LIGO noisecurve, and a = 1. 



TABLE V: Pulsars with Tdetcctabic < 3 years, assuming 100% 
conversion of spindown energy into GW energy, the Advanced 
LIGO noisecurve, and a = 1. 



Name 


Tdctoctabio/ years 


J0205+6449 

B0531+21 

J0537-6910 


2.8e+00 

5.4e-01 

2.3e+00 



energy to GW energy. For all but the millisecond pulsars 
these wave amplitudes are probably unrealistically large, 
so that this calculation gives a safe lower bound on the 
detectability timescale. We will then consider more real- 
istic scenarios, where the assumed triaxiality e is limited 
to some maximum value. 



A. Detection timescale assuming 100% conversion 
of spindown energy into GWs; probably unphysical 

We will begin by assuming the the GW and EW signals 
are perfectly in phase (i.e. a = 1 in Equation 0]). In the 
case of Initial LIGO, the (measured) timing noise and 
(estimated) GW amplitude of the Crab are both large, 
so that its timing noise detectability time is very short, 
just 0.54 years s» 200 days. However, the significantly 
lower estimated GW amplitudes of the other pulsars (see 
Figure ^) leads to only two having T^etectabie values of 
less than 3 years. The names of all three pulsars with 
^detectable < 3 years are collected in Table IIVI 

In the case of Advanced LIGO a total of 46 pulsars 
have Xdetectabie < 3 years; see Table Ivl The detectability 
timescale for the Crab is very short, just 0.058 years. 

A summary of these results is provided in the first line 
of table IVII where the number of pulsars with timing 
noise detectability timescales of less than three years are 
shown, for both Initial LIGO and Advanced LIGO. 

Setting \a\ = 10~ 2 leaves no pulsars with detectable 
timing noise using Initial LIGO, and only 7 for Advanced 
LIGO. These are listed in Table IvTTI 



B. Detection timescale assuming various upper 
bounds on the triaxiality 

The number of pulsars with timing noise detectability 
timescales less than three years, for both Initial LIGO 
and Advanced LIGO, are shown in Table IVT1 for e max = 
10~ 2 , 10~ 3 , 10" 4 , 10" 5 , 10~ 6 , 10~ 7 . For reasons of brevity 
we have not tabulated the names of the pulsars falling in 
these categories; in any case they are a subset of the 
pulsars of Tables llVl and Ivl In the case where e max takes 
on the huge value of = 10 -2 , the results are the same as 
in section lVH Al as the triaxialities required to give 100% 
GW spindown are less than this. As e max is increased, 
the number of pulsars with Tdetcctabic < 3 years decreases, 
falling to zero at e max = 10~ 7 . Clearly, timing noise will 



Name 


^detectable/ yrS 


Name 


^detectable/ yrS 


B0114+58 


2.6e+00 


J1531-5610 


2.0e+00 


J0205+6449 


1.4e-01 


B1610-50 


1.4e+00 


B0531+21 


5.8e-02 


J1617-5055 


2.8e-01 


J0537-6910 


2.3e-01 


J1637-4642 


2.2c+00 


B0540-69 


1.9e-01 


J1702-4310 


2.2c+00 


J0633+1746 


2.6e+00 


B 1706-44 


4.0e-01 


B0833-45 


2.4e-01 


B1727-33 


1.2e+00 


J0834-4159 


2.8e+00 


J1740+1000 


2.0e+00 


J0855-4644 


2.6e+00 


J1747-2958 


5.7e-01 


B0906-49 


2.7c+00 


B1757-24 


7.6e-01 


J0940-5428 


9.8e-01 


B 1800-21 


7.8e-01 


J1015-5719 


2.4e+00 


J1809-1917 


2.8e+00 


J1016-5857 


l.le+00 


B1823-13 


6.8e-01 


B1046-58 


7.3e-01 


J1828-1101 


1.9e+00 


Jl 105-6107 


9.9e-01 


B1830-08 


2.2e+00 


J1112-6103 


1.6e+00 


J1837-0604 


7.6e-01 


Jl 124-5916 


3.1e-01 


B1853+01 


2.5e+00 


B1259-63 


1.7e+00 


J1913+1011 


9.7e-01 


B1338-62 


1.6e+00 


J1930+1852 


3.0e-01 


J1420-6048 


4.3e-01 


B1930+22 


2.5e+00 


J1509-5850 


1.5e+00 


B1951+32 


5.0e-01 


B1509-58 


2.1e-01 


J2021+3651 


1.3e+00 


J1524-5625 


1.0e+00 


J2229+6114 


3.2e-01 



TABLE VI: A summary of our results for the detection 
timescale. The number of pulsars with Tdetcctabic < 3 years, 
for both Initial LIGO and Advanced LIGO with a = 1. The 
first line assumes 100% conversion of spindown energy into 
GW energy, while the remaining lines impose a cut-off in el- 
lipticity as indicated. 





Number of pulsars 


Emax 


Initial LIGO 


Advanced LIGO 


No bound 


3 


46 


10 -2 


3 


46 


10" 3 


2 


38 


10" 4 


2 


18 


10 -5 





6 


10" 6 





1 


lO" 7 









only be detectable in the GW signal if the ellipticities 
of neutron stars are 10~ 6 or larger, i.e. at least an order 
of magnitude greater than theoretical modeling currently 
suggests. 
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TABLE VII: Pulsars with Tdetectabic < 3 years, assuming 
100% conversion of spindown energy into GW energy, and 
the Advanced LIGO noisecurve, and \a\ — 10 -2 . 



Name 


Tkctcctabic/ years 


J0205+6449 


1.4e+00 


B0531+21 


5.8e-01 


J0537-6910 


2.3e+00 


B0540-69 


1.9e+00 


B0833-45 


1.7c+00 


B1509-58 


2.1e+00 


J1617-5055 


2.8e+00 



VIII. BLIND SEARCHES 

A 'blind search' is a search for gravitational waves from 
a neutron star not currently known as a pulsar. The 
presence of timing noise in the GW signal for such a 
star is more pernicious than in the case of a directed 
search. In the latter case the electromagnetic pulsar data 
can be used to demodulate the timing noise from the 
data stream. In the case of a blind search this option 
is not available, so that the timing noise will impose a 
maximum duration on any blind search beyond which the 
smooth Taylor series approximation to the phase errs by 
a significant amount (i.e. by about one radian). 

In practice, a GW astronomer would perform a blind 
search by dividing the sky up into small patches, such 
that the Doppler shifts induced by the Earth's spin and 
orbital motion are negligible over the patch |8j. Then, 
for each such patch, a search over a range of spindown 
parameters P is made. From the parametrization of 
equation (|15fl . for a given spin period P the strength 
of the timing noise is then determined, so that the deco- 
herence timescale can then be estimated from equation 
(0). The maximum search duration is then given by the 
parametrized formula: 



logT max = 2.62 + 0.93 logP - 0.531ogP_ 15 , 



(23) 



where r max is measured in years, P in seconds and P_i5 
is the dimensionless period derivative divided by 10 -15 . 
This formula could be used by GW observers to limit the 
duration of a single coherent search. 

However, there exists another limit on the length of a 
coherent blind search — namely that imposed by the re- 
quirement of doing the search in real time with finite com- 
putational power. Brady et al. [8j have investigated the 
limits placed on the length of a single non-hierarchical 
coherent search in some detail. They found that for 
a star with a GW frequency < 200 Hz and spindown 
timescale /// > 1000 years, the search is limited to a 
duration of < 18 days. Inserting these spindown param- 
eters into Equation 123(1 gives a maximum search dura- 
tion of 2.3 years, 47 times longer. For a GW frequency 
of < 1kHz and spindown timescale > 40 years, Brady et 
al find a maximum search duration of just 0.8 days; tim- 



ing noise gives a duration of 9.5 days, 12 times longer. It 
therefore appears that the finitcness of computational re- 
sources will place a more stringent limit on the duration 
of a single non-hierarchical coherent blind search than 
that due to timing noise. 

However, the high computational costs of such searches 
have motivated GW astronomers to devise more com- 
putationally efficient hierarchical search techniques. In 
these, the full data set of duration T b s is split up into 
N shorter pieces of duration T s h or t = T Q ^ S /N . Each 
short piece is then coherently analysed by matched fil- 
tering, and the separate results then combined incoher- 
ently. There are two main formulations that have been 
proposed to achieve this: the 'stack-slide' method [24| 
and the Hough transform [25|]. These will certainly be 
more computationally efficient than the single coher- 
ent analysis described above. For instance, Schutz & 
Papa 25] found that a data set of T b s = 10 7 seconds 
could be searched for GW signals of signal-to-noise ra- 
tio sa 23 by splitting the data into blocks of duration 
Tshort = 14 hours, using a 20Gflop computer, although 
this did not include a search over spin-down parameters. 

If it should prove that these hierarchical techniques are 
efficient enough that T b s should exceed the timing noise 
decoherence timescale, then the following strategy could 
be used, of a type first suggested by Brady & Creighton 
[24J in the context of variable accretion rate systems: 
Tshort could be chosen to be significantly smaller than 
T max of Equation l|23|l , so that the short coherent searches 
are not significantly degraded by the timing noise. The 
incoherent stage must then be performed allowing for 
the many possible phase shifts that timing noise could 
introduce. We suspect that this last step will be com- 
putationally expensive as the timing noise can map out 
many different paths through phase space, but we won't 
attempt to quantify this cost here. 



IX. CONCLUSIONS 

In this paper we have addressed three issues. 

First we asked for which of the known pulsars might 
timing noise prevent a GW detection, if the experimenter 
assumes only smooth Taylor-series spindown. Our main 
conclusions were as follows: 

• If the GW timing noise is at the same level as the 
EW timing noise (i.e. a = 1 in Eqn. 0]), then, 
assuming 100% conversion of spindown energy into 
GW energy, as many as 3 pulsars may be rendered 
undetectable by Initial LIGO (Table [n|, and as 
many as 9 for Advanced LIGO ( Tabic liTTjl . 

• If the GW timing noise is weaker, at the level 
of |cv| = 10~ 2 , no pulsars will be rendered unde- 
tectable by timing noise. 

• Dropping the assumption of 100% conversion of 
spindown energy into GW energy and instead plac- 
ing an upper bound on the ellipticities leads to 



fewer stars for which decoherence may be an is- 
sue. Setting e max = 1CP 7 or smaller we find that 
timing noise will not prevent the detection of any 
pulsars. 

Next we asked for which of the known pulsars might 
timing noise be strong enough to be observed in the GW 
data stream. Our main conclusions were as follows: 

• If the GW timing noise is at the same level as 
the EW timing noise (i.e. a — 1), then, assum- 
ing 100% conversion of spindown energy into GW 
energy, timing noise will be detectable in 3 pulsars 
by Initial LIGO (Table HVj), and in 46 pulsars for 
Advanced LIGO (Table EJ. 

• If the GW timing noise is weaker, at the level of 
|a| = 10~ 2 , timing noise will not be detectable in 
the GW signal of any pulsars using Initial LIGO, 
and in only 7 using Advanced LIGO (Table |V1*T|) . 

• Dropping the assumption of 100% conversion of 
spindown energy into GW energy and instead plac- 
ing an upper bound on the ellipticities leads to 
fewer stars for which timing noise is detectable (Ta- 
ble |JT}. Setting e max = 10~ 7 or smaller we find 
that timing noise is not detectable in the GW sig- 
nal of any pulsars. 

Finally we asked how timing noise might affect blind 
GW searches, i.e. searches for GW emitters not currently 
observed as pulsars. Our main conclusion was that: 

• Timing noise places a limit on the length of a data 
set that can be coherently analyzed f Equation I23[l . 



Proposed hierarchical search techniques could per- 
form searches over longer durations by allowing for 
the possibility of (many different possible realisa- 
tions of) timing noise when the short coherent anal- 
yses are incoherently combined. 



Table [H] is probably of most interest to today's GW 
data analysts — it lists the pulsars for which timing noise 
is most important for the first generation interferometers. 
For these three pulsars, GW astronomers may wish to use 
electromagnetic timing residuals to improve their ability 
to search for these stars. A method of doing so for one 
of these (the Crab) was described recently by Pitkin & 
Woan 0. 

To sum up, timing noise may be an important feature 
of the GW signal of some tens of young pulsars, but these 
stars must have very large ellipticities in order for the 
GW emission to be strong enough to be detectable. If 
observed in the GW signal, timing noise would provide a 
new insight into neutron star dynamics. 
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